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We outline the basic properties of a pertubative QCD factorization formalism that maintains exact 
over-all kinematics in both the initial and final states. Such a treatment requires the use of non- 
perturbative factors that depend on all components of parton four-momentum. These objects 
are referred to as parton correlation functions. We describe the complications faced in defining 
parton correlation functions and discuss recent progress. Emphasis is placed on the need for 
precise operator definitions in a complete derivation of factorization. 
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1. Exact over-all kinematics and the need for a generalized treatment of 
factorization 

The standard collinear factorization theorems of perturbative QCD (pQCD) rely on a number 
of kinematical approximations that change the momentum of the final state particles. In inclusive 
deep inelastic scattering (DIS), for example, the struck parton is assumed to have a non-zero com- 
ponent of four-momentum only in the plus direction corresponding to the target beam direction. 

It has been known for some time that neglecting the transverse component of the struck par- 
ton momentum is not always valid. This has lead to studies of "^^ -factorization" involving kf- 
unintegrated (or just unintegrated) parton distribution functions [^. More recently, it has been 
noted that in certain cases over-all four-momentum conservation (involving both transverse mo- 
mentum and paiton invariant energy) must be enforced to avoid making large errors, particularly 
when the details of final states are important [Q]. This has motivated the formulation of factor- 
ization theorems in which no approximation is made on the momenta of initial and final states. 
The non-perturbative objects in such a factorization formula will depend on all components of par- 
ton four-momentum. We call these fully unintegrated objects parton correlation functions (PCFs). 
In these proceedings we outline the basic structure of the fully unintegrated formalism proposed 
in [^]. In addition, we re-emphasize the need for exact operator definitions for the non-perturbative 
factors in a factorization formalism. 

2. The role of operator-based definitions 

Having non-perturbative factors rooted in operator definitions is important for the derivation 
of a reliable factorization formula. To understand this, let us briefly review the basic requirements 
of a factorization theorem: 

• For a given process with hard scale Q, a factorization formula exists if the cross section can 
be written approximately as a generalized product of several factors. The hard scattering 
coefficient C should involve only lines that are off-shell by order Q and can be calculated 
explicitly in pQCD. The other factors involve infrared and collinear lines and parameterize 
the non-perturbative physics. Errors should be suppressed by powers of A/2 where A is a 
typical hadronic mass scale. 

• The non-perturbative factors must be parameterized by experimental data. But if factoriza- 
tion formulae exist for different processes, and involve the same non-perturbative factors, 
then one can parameterize a non-perturbative factor in one experiment and use it in another 
to make first-principle predictions. If this can be shown, we say that the non-perturbative 
factors are universal. 

Operator definitions are what allow for a comparison of the soft and collinear factors used in differ- 
ent processes. Hence, they are needed if one is to have confidence in the second buUeted statement 
above. As an example, consider the simple case of totally inclusive DIS. The definition of the 
integrated parton distribution function should arise naturally from the sequence of approximations 
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needed to factorize the hadronic tensor. The hadronic tensor written using standard notation is, 

W^\P,q) =4K'l^{p\J^mX){X\rmp)5^'\px - ip + q)). (2.1) 

X 

Starting with this very basic formula and applying certain approximations, one can derive the for- 
mula, 

P^yW^" {P, q) = P^.H^" (<7, § , At ) ® /y/, , M ) + ^(A/e IPuvW^" (P, q)\), (2.2) 

where P^v is defined to act on the electromagnetic vertices to project out a particular structure 
function, say Pi. The function /y/;,((§,At) is the usual expression for the fully integrated parton 
distribution function (PDF), 

fjd^,^) = I X {p\m^-,OT)vUuj)r^Voiu,Mo)\p)- (2.3) 

Here \j/ is the quark field operator, and Vq («/) is a Wilson line operator in the direction = (0, 1 , 0/ ) 
needed to make the definition exactly gauge invariant. The lowest order hard scattering matrix 



element is the usual one involving only the electromagnetic vertex. Using it in Eq. ( |2.2| ) reproduces 
the parton model. Higher order corrections are calculated by considering more complex graphs and 
applying a sequence of subtractions to remove double counting. The unregulated PDFs contains 
the usual ultraviolet divergences which are effectively removed using standard renormalization 
techniques, with a renormalization scale /i. The resulting evolution equations describe the well- 
known scaling violations of DIS. The fact that the PDF ,At) also appears in the factorization 
formula for other processes means that it can be parameterized and used for future predictions. 

One may hope to have the same powerful structure in a more general unintegrated formalism. 
In a treatment that includes the transverse momentum of the struck parton, it is tempting to propose 
a differential definition for the -unintegrated parton distribution. 



(2.4) 



Studies of small-x physics ^ suggest that this is probably correct to a good approximation in the 
small-A; limit. However, it is unclear in general whether there is a reliable sequence of approx- 



imations that allow ( |2.4[ ) to be factorized out of the hadronic tensor with only a hard scattering 
coefficient left over. The universality of the -dependent PDF is thus called into question. 

The situation with /c, -unintegrated PDFs is further complicated by the fact that the most obvi- 
ous candidate for a definition is unsuitable for use in factorization. Namely, if we leave one of the 



integrals in Eq. (2.3) undone, we obtain the seemingly natural definition. 



^j/p{x,kT,^) = I ^ ^p\^j(^o,w-,WT)V:{n)r+Vo{n)YjmP)- (2-5) 

However, this definition acquires divergences from gluons moving with infinite rapidity in the 
outgoing quark direction making it inappropriate for use as a PDF. (See [Q] for a more detailed 
discussion.) An additional problem, pointed out by Belitsky et al. |^, is that exact gauge invariance 
requires the insertion of a gauge link at light-cone infinity. Recent work in providing consistent 
operator definitions for ^, -unintegrated PDFs has been done by Hautmann and Soper [0]. 
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(a) (b) 
Figure 1: (a) General contribution to DIS at lowest order and (b) Structure of the graph after factorization. 

3. A fully unintegrated formalism 

As already mentioned, there are situations where we need a formalism based on non-perturbative 
objects that depend on all components of parton four-momentum. An approach to factorization us- 
ing PCFs was formulated in [js]] for the case of a scalar field theory, and extended to the case of 
a gauge theory in [^.^. We summarize these results now. To save space we do not list the actual 
definitions, but rather outline the basic structure of their derivation. For all details, the reader is 
referred to |^. 

To avoid making errors of the type discussed in the first section, one must begin with graphs 
of the general structure shown in Fig. |l|(a) rather than the usual handbag diagram. The bubbles 
represent sums of diagrams contributing to initial and final states. The extra gluons shown attaching 
the colUnear bubbles to the hard scattering bubble represent possible extra target and final state 
collinear gluons. In addition, there may be arbitrarily many soft gluons connecting the outgoing 
jets via a soft bubble. Before factorization, these bubbles are <I> for the target collinear lines, ^ 
for the jet collinear lines, and ^ for the soft lines. 

It follows from general arguments [0] that graphs with the topology of Fig. [I](a) are the most 
general contributions to DIS involving a single outgoing jet. Applying Ward identities to the sum 
of graphs with this allows the extra soft and collinear gluons to be disentangled into separate fac- 
tors. After topological factorization is achieved, the graph takes the form shown in Fig. |l|(b). The 
different PCFs include a soft factor S, a jet factor J, and a target factor(fully unintegrated PDF) F. 
The PCFs are represented graphically by the bubbles. The double lines associated with each bubble 
are eikonal lines that correspond to Wilson lines in the definitions of the PCFs. To be consistent 
with factorization, the definitions of the PCFs also require double counting subtractions. (To avoid 
clutter, the double counting subtractions aren't shown explicitly in the figure.) The factors S,J, and 
F can ultimately be shown to arise from operator definitions of the PCFs. Schematically, the final 
factorization formula is, 

a = C(g)F®7®5+^((A/er|a|), (3.1) 
where (7 is a measurable quantity such as a cross section or structure function and a > 0. 

'The derivation is so far only complete for the case of an abelian gauge theory 
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4. Outlook 



Our factorization formula is now more complicated than the usual one in Eq. (2.2). Each 
PCF depends on several parameters and each needs to be fitted to experimental data. For the fully 
unintegrated formalism to be practical, factorization formulae using the same PCFs will need to be 
derived for a number of non-trivial processes. In addition to evolution in /i, the evolution equations 
for the PCFs will also involve evolution in other rapidity variables which act as effective cutoffs on 
rapidity divergences. One hope is to relate this type of evolution to more common approaches such 
as the CCFM equation. 

As we have mentioned, the factorization formula represented schematically in Fig. [l](b) is only 
complete at lowest order in the hard scattering coefficient because it involves only one outgoing jet 
line. However, this result is already quite useful because it now allows higher order corrections to 
be obtained via double counting subtractions applied to more complicated graphs. As in for the 
scalar field theory, the hard scattering coefficients are expected to be ordinary functions as opposed 
to the generalized functions (e.g. 5-functions) that appear in the hard coefficients of the standard 
integrated formalism. 
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